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Abstract: Let X = {Xt)t>o be a Levy process with absolutely continu- 
ous Levy measure u. Small time polynomial expansions of order n in t are 
obtained for the tails P {Xt > y) of the process, assuming smoothness con- 
ditions on the Levy density away from the origin. By imposing additional 
regularity conditions on the transition density pt of Xt , an explicit expres- 
sion for the remainder of the approximation is also given. As a byproduct, 
polynomial expansions of order n in t are derived for the transition densities 
of the process. The conditions imposed on pt require that its derivatives 
remain uniformly bounded away from the origin, as t — > 0; such conditions 
are shown to be satisfied for symmetric stable Levy processes as well as for 
other related Levy processes of relevance in mathematical finance. The ex- 
pansions seem to correct asymptotics previously reported in the literature. 
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1. Introduction 

Levy processes are important building blocks in stochastic models whose evo- 
lution in time might exhibit sudden changes in value. Such models can be con- 
structed in rather general ways, such as stochastic differential equations driven 
by Levy processes or time-changes of Levy processes. Many of these models have 
been suggested and heavily studied in the area of mathematical finance (see [3] 
for an introduction to some of these applications) . 

A Levy process X — {Xt)t>o is typically described in terms of a triplet 
(cr^ , b, v) such that the process can be understood as the superposition of a 
Brownian motion with drift, say aWt + bt, and a pure-jump component, whose 
discontinuities are determined by v in that, the average intensity (per unit time) 
of jumps whose size fall in a given set of values A is v{A). Thus, for instance, if 
^"((—00, 0]) = 0, then X will exhibit only positive jumps. A common assumption 
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an interesting counterexample. 
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in many applications is that ly is determined by a function s : M\{0} — > [0, c»), 
called the Levy density, in the sense that 

„{A) / s{x)dx, yA e 6(M\{0}). 



J A 

Intuitively, the value of s at xq provides information on the frequency of jumps 
with sizes "close" to xq. 

Still, Levy models have some important shortcomings for certain applications. 
For instance, given that typically the law of Xt is specified via its characteristic 
function 

neither its density function pt nor its distribution function P (Xt < y) are explic- 
itly given in many cases. Therefore, the computation of such quantities neces- 
sitates numerical or analytical approximation methods. In this paper we study, 
short time, analytical approximations for the tail distributions P {Xt > y). This 
type of asymptotic results plays an important role in the non-parametric esti- 
mation of the Levy measure based on high-frequency sampling observations of 
the process as carefully reported in [5] (see also [T7], H], and [ID])- In Section [2j 
we present some of the ideas behind this important application of our results. 

It is a well-known fact that the first order approximation is given by tv(\y, oo)), 
in the sense that 

f^^-^V{Xt>v) = v{[y,^)), (1.1) 
provided that y is & point of continuity of v (see, e.g.. Chapter 1 of Bertoin [1]). 



A natural question is then to determine the rate of convergence in ( 1.1 1. In case 
of a compound Poisson process, this rate is 0{t), and it is then natural to ask 
whether or not the limit below exists for general Levy processes: 

lim J|^P(Xt>y)-K[y,oo))|. (1.2) 

In this paper, we study the validity of the more general polynomial expansion: 

" 4-k 4-n+l 

V{Xt>y) = y^dk-^+—nr,{t), (1.3) 
fe=i 

for certain constants dk and a remainder term TZn{t) bounded for t small enough. 



Note that in terms of the coefficients of (1.3), the limit (1.2 1 is given by 



^ = lim \ {\v{Xt >y)- iy{[y,^)) } . (1.4) 



For a compound Poisson process, the expansion ( 1.3 1 results easily from con- 
ditioning on the number of jumps on [0, t\. Thus, infinite-jump activity processes 
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are the interesting cases. Ruschcndorf and Woerner [TS] (see Theorem 2 in Sec- 
tion 3) report that for a fixed > 1 and rj > 0, there exists a s'{N) > and 
to > such that, for all e G (0, e'{N)) and t e (0, to), 



N-l 

P(^t>y)= V ^i^r([2/,oo)) + 0,,^(t^), fory>7^, (1.5) 



i=l 



where v^{dx) = l{|^|>£j.z^((ia;). When A'^ = 3, this result would imply that, for 
< e < ?;/2Ae'(iV),~ 



> y) ^tv{[y,^)) + ^ I I l{„+„>y}j.(rff)i/(du) + 0,,„{t') 



\u\>€ J \v\>e 



Thus, (1.5) would imply that 

]^'^^\\\'^{Xt>y)-v{[y,co))\ = \( I l{,,+.^>y}v{dv)v{du), 

I I r J ^ J\iL\>e J\v\>e 

which is independent of the Brownian component aWt and of the "drift" bt. We 
actually found that this limiting value is not the correct one and provide below 
the correction using two different approaches. Let us point out where we believe 
the arguments of |18j are lacking. The main problem arises from the application 
of their Lemma 3 in Theorem 2 (see also Lemma 1 in Theorem 1). In those 
lemmas, the value of actually depends on 5. Later on in the proofs, 5 is taken 
arbitrarily small, which is likely to result in to ~^ (unless otherwise proved). 

We prove ( 1.3 1 using two approaches. The first approach is similar in spirit to 
that in [TH]. It consists in decomposing the Levy process X into two processes, 
one compound Poisson process X'^ collecting the "big" jumps and another pro- 
cess Xf accounting for the "small" jumps. By conditioning on the number of 
big jumps during the time interval [0,t], it yields an expression of the form 

P {X, >y) = e-'-' (xt + J2^^>y]+ Oir+'). 

k=l ■ V j=l / 

By taking a compound Poisson process X^ with jumps {£,i\i>i having a smooth 
density, one can expand further each term on the right-hand side using the 
following power series expansion: 

^ j-k +n+l /•! 

E5(X,) = 5(0)+^^i'=5(0) + — / (l-a)"E {L"+ig(X„,)}d«, (1.6) 
fe=i ■ -^^ 

valid for any n > and g £ C^"^^, the class of functions having continuous 
and bounded derivatives of order Q < k <2n + 2. Above, L is the infinitesimal 
generator of the Levy process, i.e., 

{Lg){x) ■■= ^g"{x) + bg'{x)+ / {g{u + x) - g{x) - ug' {x)l{\u\<i}) v{du), 

(1.7) 
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for any function g G 



For n — 0, (1.6) takes a familiar form (see e.g. Lemma 19.21 in [H]): 



Eg{Xt)^g{0) + t E {Lg{X^t)}da = g{0) + E{Lg{Xu)}du, (1.8) 
Jo Jo 

which is an easy consequence of I to's formula. The general cas e fol lows eas- 
ily by induction in n. Indeed, if (1.6 1 is valid for n, applying (1.8 1, /q^(1 — 
a)"E {U'+^g{Xat)} da becomes 



(1 - a)" <^ L"+\g(0) + at 



1 



n+1 



rn+l 



5(0) 



n+1 



E {L"+^g{Xa,,o.t)}da'^da 
(1 - d)"+i E {L''+'^g{X&t)} da, 



where we changed variables a := aa' and applied Fubini' s Theorem. Another 



proof of (1.6 1 is given in [9] based on Fourier approximations of g (Proposition 
1 and 4 in there). 

In the second order approximation case (rt = 2), we give another proof for 



(1.3 1 which relaxes the assumptions on the Levy density s, by requiring only 



smoothness in a neighborhood of y and local boundedness away from the origin. 
This approach is based on the following recent asymptotic result by Jacod [TO] : 



linii Eg{Xt)=<T'^ + j g{x)v{dx), 



(1.9) 



valid for a j/-continuous bounded function g such that g{x) ^ x^, as a; — > 0. 
In case the process is of finite variation and has no diffusion term, we prove 
the second order expansion as long as s is continuous at y and locally bounded 
away from 0. We also present a counterexample, originally suggested by Philippe 
Marchal, which shows that the result is not valid if s is not continuous (see |13j 
for further developments). 



In order to provide explicit formulas for the coefficients dk in (1.3 1, in Section 



|4] we consider a second approach whose basic first step is to approximate the 
indicator function l[j,,oo) by smooth functions /„j in such a way that 

lim Ef„,{Xt)=V{Xt>y). 

m — ^oo 



The idea to derive ( 1.3 ) is to apply ( 1.6 1 to each smooth approximation /, 



and 

■ 00. 



show that the limit of each term in the power expansion converges as m 
We emphasize that this approach is carried out without additional assumption 
on s, except smoothness and local boundedness away from the origin. In Section 
[5] we exploit further the approximation of / by the smooth functions fm to 
provide an explicit formula for the remainder TZn{t) in (1.3 1. To carry out this 
plan, we impose more stringent conditions on X than those required in the 
first approach. In particular, we require that Xt has a C°°-transition density pt, 
whose derivatives remain uniformly bounded away from the origin, as i — + 0. 
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As a byproduct of the explicit remainder, polynomial expansions of order n in t 
are derived for the transition densities of the process extending a result in |TH]. 

In Section [6] the boundedness conditions on the derivatives of the transition 
densities are shown to hold for symmetric stable Levy processes. The validity of 
this uniform boundedness for general tempered stable processes is also considered 
in Section [7] via a recursive formula for the derivatives of the transition density. 
Tempered stable processes have received a great dealt of attention in the last 
decade due to their applications in mathematical finance. Among their members, 
we can list the CGMY model of |2] . See Rosihski for a detailed study of this 
class of processes. 

We note finally, that throughout the paper we only consider asymptotics 
for F{Xt > y),y > 0, but that our methodology also gives results for P{Xt < 
-y),y > 0, replacing j/([y,+oo)) by j/((-oo, -j/]). 



2. An application: nonparametric estimation of the Levy density 



In this part we present an application of the small-time asymptotics considered 
in this work as a matter of motivation. One problem that has received attention 
in recent years is that of estimating the Levy density s of the process in a non- 
parametric fashion. This means that, by only imposing qualitative constraints on 
the Levy density (e.g. smoothness, monotonicity, etc.), we aim at constructing 
a function s that is consistent with the available observations of the process X. 
The minimal desirable requirement of our estimator s is consistency; namely, 
the convergence s — > s, say in a mean-square error sense, must be ensured when 
the available sample of the process increases. 

When the data available consists of the whole trajectory of the process dur- 
ing a time interval [0,T], the problem is equivalent to estimating the intensity 
function of an inhomogeneous Poisson process (see e.g. [TS] for the case of finite 
intensity functions and for the case of Levy processes, where the intensity 
function could be infinite). However, a continuous-time sampling is not feasi- 
ble in reality, and thus, the relevant problem is that of estimating s based on 
discrete sample data Xt^^, . . . ,Xt^ during a time interval [0,T]. In that case, 
the jumps are latent variables whose statistical properties can in principle be 
assessed if the frequency and time horizon of observations increase to infinity. 



It turns out that asymptotic results such as (1.2 1 and (1.3 1 play important 
roles in determining how frequently one should sample (given the time horizon T 
at hand) such that the resulting discrete sample contains sufficient information 
about the whole path. We can say that a given discrete sample scheme is good 
enough if we can devise a discrete-based estimator for the parameter of interest 
that enjoys a rate of convergence comparable to that of a good continuous- 
based estimator. Let us explain this point with a concrete example. Consider 
the estimation of the following functional of s: 

(3{(p) :— / (p{x)s{x)dx, 
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where iy9 is a function that is smooth on its support. Assume also that the support 
of if is an interval [c,d] so that the indicator l[c^d] vanishes in a neighborhood 
of the origin. A natural continuous-based estimator of (3{ip) is given by 

s<T 

Using the well-known formulas for the mean and variance of Poisson integrals 
(see e.g. [19, Proposition 19.5]), the above estimator can be seen to converge to 
P{(p), and moreover, 

We can thus say that P^^if) converges to f3{ip) at the rate of 0(T~^/^), in the 
mean-square sense. 

Suppose that instead we use a reasonable discrete-based proxy of /3^, using 
the increments Xt-^ — Xtg, . . . , Xt^ — Xt^_-^ of the process instead of the jumps 

1 " 

i=l 

where ir : to < ■ ■ ■ < tn = T . A natural question is then the following: How 
frequently should the process be sampled so that (3'!^{(p) (3{<f) at a rate of 
0(T-i/2)? To show in a simple manner the connection between the previous 
question and the asymptotics (1.2 1, suppose that the sampling is "regular" in 
time with fixed time span A„ := T/n between consecutive observations. In that 
case, we have 

E (/3J(^) - (3{^)r < i/3(v.2) + 1 I i- E ^2 (^^^^ ) _ ^(^2) 

1 

— E^(XaJ-/3(^) 
From the previous inequality we see that the rate of convergence in the limit 

^m^E^{XA)^f3{v), (2.1) 

will determine the rate of convergence of P^iffi) towards f3{(p). To determine the 
rate of convergence in ( |2.1| , one can simply link E(/j(Xa) to P(Xa > y), and 
link P{ip) to iy([y, oo)). This is easy if ip is smooth on its support [c,d]. Indeed, 
we have that 



A 



< (ll<y=lloo + |l<^'|ll) sup 
yelc,d\ 



^P[^A>2/]-K[2/,oo)) 



Hence, the rate of convergence of A ^¥{Xa > y) towards I'dy, oo)) determines 
the rate of convergence of A~^ E ip {X^) towards P{ip). In particular, the result 
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(O) will tell us that, for PJ^{(p) to converge to P{ip) at a rate of OiT-^/"^), 
in the mean-square sense, it suffices that the time span b etwe en consecutive 
observations A is o(r~^/^). ft is important to remark that (1.2 1 can be seen to 
hold uniformly in y > y, for an arbitrary y > 0. 



The ideas outlined in this section, as well as the asymptotic result (1.2 1, are 
heavily exploited in |3] and [5], where the general problem of nonparametric 
estimation of the Levy density s is studied using Grenander's method of sieves. 



3. Expansions for the transition distribution 

As often, e.g. see [TB], the general strategy is to decompose the Levy process into 
two processes: one accounting for the "small" jumps and a compound Poisson 
process collecting the "big" jumps. Concretely, suppose that X has Levy triplet 
(cr^,&, j^); that is, X admits the decomposition 

Xt = bt + aWt+ I I x{n- p,){dx,ds) + [ [ Xfi{dx,ds), (3.1) 

^0 J\x\<l Jo J\x\>l 

where 14^ is a standard Brownian motion and fi is an independent Poisson mea- 
sure on M_|_ X M\{0} with mean measure fl{dx, dt) := v{dx)dt. Note that is 
the random measure associated to the jumps of X. Given a smooth truncation 
function Cg € C°° such that l[e/2,e/2](a;) < c^{x) < l[^^^]{x), set 



X^ := / / xce{x)^i{dx,ds), (3.2) 
Jo Jr 

XI X, - Xt, (3.3) 

where Ce{x) := 1 — Cg. It is well-known that X'^ is a compound Poisson process 
with intensity of jumps := J Cj(x)z^(dx), and jumps distribution Cj(a;)z^((ia;)/Ae 
The remaining process X^ is then a Levy process with jumps bounded by e and 
Levy triplet (cr^, &£, Ce(a:)i^((ix)), where 

b^:=b— / xc^{x)v{dx). 

J\x\<l 



There are two key results that will be used to arrive to ( 1.3 1. The first is the 



expansion (1.6 1. The following tail estimate will also play an important role in 
the sequel: 

V{\XI\ >y)< exp{a2/ologyo}exp{ay-aylogy}ty°, (3.4) 

valid for an arbitrary, but fixed, positive real a in (0,e~^), and for any t,y > 
such that t < i/Q^y, where yo depends only upon a (see [THl Lemma 3.2] or [ini 
Section 26] for a proof) . 
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Remark 3.1. For an alternative proof of (3.4-), use a generic concentration 
inequality such as ^8, Corollary 1] to get (when a = Q): 

nxt>y)^nxt-¥,xi>x) 

whenever x := y — E^f > 0, and with :— J^^^^^u'^i'{du) . Now EXf = 

t{b^ + Xri<|3.|<£} xv{dx)), and ast Q, x ^ y and {eV'^ /exY^^ — s- {eV"^ / eyY ^ ^ , 

with moreover t^l^ jt^ = cxp{{y — EXf — 2e)logt/e) 0, as long as y > 2e. 
Finally, since as t — > 0, P((7VKt > y/2)/t'^ 0, the general case follows. 



We are ready to show (1.3 1. Below, is the infinitesimal generator of X'^ 
and we use the following notation: 

Se := CeS, Se:=l-Se, L"g^g, s*^ = 

s»{x) = I -sf'^\x - u)Uu)du, {i > 2), sf*g^ g. 



Theorem 3.2. Let y > 0, n > 1, and < e < y/{n + 1) A 1. Assume that v 
has a density s such that for any < k < 2n + I and any S > 0, 

Qk.s ■= sup \s^''Hx)\ < oo. 

\x\>S 

Then, there exists a to > such that, for any y > y and < t < to, 

" fj 

V{Xt > y) = e-^^*^c, - + Oe,,(r+i), (3.5) 

3 = 1 



J , 

J 



where 

i=l ^ 

with fi{x) := J^^s*/{u)du. 

Proof. Throughout this part, we write f{x) := l[x>y}- In terms of the decompo- 
sition X := X'^ + X"^ described at the beginning of this section, by conditioning 
on the number of jumps of X^ during the interval [0,t], we have that 



oo 



E/(XO- E/(Xne-^^* + e-^^* ^ X.^ + ^C. (3.6) 



k=n+l 



-"'^*E^IE/ + (3.7) 

fc=i ■ V 1=1 / 
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where Ce{x)s(x)dx/ X^. Taking a :— (n + l)/y, (3.4) and < f < I imply 
that the two terms on the right hand side of (3.6 1 are Oe,y{t"'^^) as t — > 0, 
provided that t < to -.^ i/Q^y. Next, for each fc > 1, 



where 



fk{x) E/ x + ^e. = 



e=i 



a; + ^e. > y 

\ 1=1 ) 



which is Cf'+^ since the density of is C^"^^. Then, one can apply to 
get 



^ ^ y-i j.n-\-\ — k /*! 



1=0 



(3.8) 



Let Le be the infinitesimal generator of X'^ , given by 

{L,g){x)^h,g\x) + ''^g"{x)+ J ^ g" {x + Pw){l - /3)d(3w^c,{w)s{w)dw, 

for g e C^, and for k> 1, let 

dnf. n^=i(l - (3i)dPew]ce{wi)s{we)dwi, 
which clearly a finite measure on [0, 1]'' X M*^. Then, note that 



{Llg){x)^ J2 cJl)Alg{. 



(3.9) 



where /Cj :— {k :— (fci, ^2,^:3) : ki + k2 + — i}, 



fe3 



(fcl+2fc2+2fc3) 



if fcs > 1 and Akg(x) :== 5('=i+2'=2)(a:), if /sa = 0. Since 

fi'\x)^X^''i~ir'-s<''-'K-s(t'Hy~x), 

and Se(-) G C^"^^, there exists a constant &„^e < 00 (independent of y), such 



that 



/fell (a2n+l,e/2) 
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and so, the last term in (3.8 1 is 0(i"+^ ''). Plugging (3.8) into (3.6 1 and rear- 
ranging terms, we get 



n , ^ m / 



j=l ^- k=l 



which is exactly (3^), because A^'i™ ''ft = L'^ ^ fk- 
Remark 3.3. 



□ 



(i) The expansion (1.3) follows from (3.5). Indeed, expanding e we get 
that for any y > y and < t < to: 



k=l 



nth 



j=i v-'/ 

In the next section we give a more explicit expression for dk. 
(a) The first two terms in { 3.1lJj^ can be easily computed: 



(3.10) 



(3.11) 



di — s{u)du ^ h'{[y,(X))) 
■'y 

^2 = -2AsI^([y, OO)) + // l{ui+u2>y}Se{ui)Se{u2)dUidU2 



a^s'{y) + 2b^s{y) ^2 / s'{y ~ /3w){l - (3)dPw^Se{w)dw. 



(Hi) The coefficients dk in {3.10) are independent of e since they can he defined 
iteratively as limits ofV{Xt > y). For instance, 

lini \v{Xt >y)^ di, lim \ \-V{Xt > y) - d\ = da- 

One can obtain an expression for d2 that is independent of e by taking 
the limit as e — )■ 0. For instance, if X is of bounded variation with drift 
bo :— b — J|^|<;]^ xv{dx) and volatility a, then ^2 becomes 



d2 = -a's'iy) + 2bos{y) - {v{[y, <x^))Y 



V ry 



s(u)dus{x)dx + 2 



J y—x 



oo fO 



s{u)dus{x)dx. 



y ''y-x 
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In general, it turns out (see the Appendix) that 6,2 "simplifies" to the 
following expression when e — )■ 0; 

d2 = -o^ s\y) + 2hs{v) - v(iy, + i.((2//2, y)f 

/ s{u)dus{x)dx — 2s(y) / xs{x)dx 

-oc Jy~x Jy/2<\x\<l 

/ {s{u) - s{y)}dus{x)dx. □ 

-y/2 Jy-x 



We now present an alternative proof for the expansion (1.4) that requires less 
stringent assumptions. The following asymptotic result due to Jacod [10^ will 
be of importance: 



(3.12) 



lim- Eg{Xt)^a^+ / g{x)iy{dx), 
t— >o t J 

valid if g is z^-continuous, bounded, and such that g{x) ~ x^, as a; — *■ 0. 

Proposition 3.4. Let y > and < e < y/2Al. Assume that v has a density s 
which is bounded outside of the interval [— e,e], and that is in a neighborhood 
ofy. Then, the limit (H) exists and can be written as: 

do / , ^ . , ^ 

f = -y^'(y) + &e.<2/) + 



{s{y - u) - s(?/)} duse{x)dx 

'^{x+u>y}Ss{u)Se{x)dudx ~ Xel^{[y, Oo)). 



Proof. Let f{x) := l^^yyy and let 



Ait) := 



1 fl 



t [t 

With the notation of Theorem 13.21 we have 



EfiXt)- / fixHdx) 



A(<)= iE/(xne- 



1 



{Ef{X^ + x)- f{x)}s,{x)dx (3.13) 



1 -e 



f(x)Ux)dx-re~^^'Y. -^^ E/ + , 



n=2 



since e < y/2. In view of (3.4|, the first term on the right hand side vanishes 
when i — > 0. Then, except for the second term, all the other terms are easily 
seen to be convergent. Let us thus analyze the second term. Let 



{Ef{Xl + x)- f{x)}-s,{x)dx. 
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Since < e < y/2 and the support of q is [— e,£], B{t) can be decomposed as 

ry rv+s 
B{t) := P {XI >y-x) s{x)dx - / P {X^ < y - x) s{x)dx 



P {X^ >y~x) Se{x)dx - P{XI <y-x} Se{x)dx. 

x<y—e J y+e 



Since s is bounded and integrable away from the origin, the last two terms can 
be upper bounded by A^P {\Xl | > e} , which, divided by t, converges to in view 
of ( 1.1 1. After changing variables to u = y — x and applying Fubini's Theorem, 



the first term above becomes: 
ry 



(a;Ae)VO 



{XI >y~x) s{x)dx = / P {XI > u) s{y - u)du = E /+ {XI) , 



where f+{x) :— s{y — u)du. Similarly, 



' {XI <y-x) s{x)dx = / P {XI < -u) s{y + u)du = E /_ {XI) , 

"'0 

where f-{x) := ^^^^"^^ s{y + u)du. Next, consider the function 

f+{x) - s{y) {xAe), x>0 

f{x) : 

-f-{x) + s{y){-x he), x<0, 



and note that lima;^o f{x)/x'^ = —s'{y)/2. In view of (3.121, we conclude that 

lim^E/(Xf) = -^a2 + J f{x)c,{x)s{x)dx. 

Thus, the sum of the first two terms in the decomposition of B{t) are 

¥.f+{XI)-¥.f_{XI)^¥.J{XI) + s{y)¥,h{XI), (3.14) 



where h{x) = a;l|j.|<g — e\x<-e + ela;>e- Let us analyze the last term in (3.141: 

lim \v.h{XI) ^ lim \kXI~ Imi ^ EXf l{|xf |>e} 

+ e lim -P{XI > e} - e lim -P{Xf < -e) h^. 
We are finally able to give the limit of B{t)/t: 

Imi -^B{t) ^ - ^ ^2 + s{y) b, + I f{x)c,{x)s{x)dx. 
A little extra work leads to the expression in the statement of the result. □ 
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It is not clear whether or not Proposition |3.4| remains true when a ^ and 
the density of ly is not differential in a neighborhood of y.li a = 0, one can relax 
the differentiability condition as follows. 

Proposition 3.5. Let y > and < e < y/2 A 1. Assume that v has a density 
s which is bounded outside of the interval [— e, s] and that is continuous in a 
neighborhood of y. Assume also that cr = and that 



{bl<i} 



|a;|i/(c?a;) < oo. 



(3.15) 



Then, the limit (1.2) exists and is given by 

s(y — u)duSi;{x)dx 



y := bo s{y) + 



'^{x+u>y}Se{u)Se{x)dudx - \eV{[y, Oo)), 



where bo := b — Ji^ki xv{dx) 



Proof. The proof is very similar to that of Proposition 3.4 However, instead of 



(3.121, we use the following asymptotic result 

lim^ Eg(Xt) = |6o| + j g{x)v{dx), 



(3.16) 



valid for any continuous bounded function g such that g{x) ~ | (see 

e.g. Jacod [1^). Define the function 



f+{x), x>0, 
-.f-{x), x<0, 



and note that lima;_fO f{x)/x = s{y). By (3.16), 



lim-E/(Xf) = s(y)6o + 



f{x)v{dx). 



(3.17) 



|2:|<e 



Using the arguments of the proof of Proposition 3.4 (3.17) implies that 



}^l\B{t)^bos{y) + 



s{y - u)ss{x)dx. 



and this gives the value of c?2 stated in the statement of the result. □ 

Remark 3.6. The continuity of s is needed in Proposition as the fol- 
lowing example suggested by Philippe Marchal shows (see llSJ, for further de- 
velopments). Let Xt :— St + Yt, where S is a strictly a-stable Levy process 

such that St^t^^°'Si and Y is an independent compound Poisson with jumps 
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and jump intensity 1. Suppose that 1 < a < 2 and that the density p of 
^ is such that p{y~) ^ pijj'^). Notice that the Levy density of the process is 
s{x) — x~°'~^ +p{x). Then, as shown next, 

C{t) := ^ |ip(X, >y)- ^.([y,oo))| , 



converges to a non-zero limit as t —> 0, and so, (1.2) is infinite. Indeed, condi 



tioning in the number of jumps of the compound Poisson component Y^, 

oo „ / n \ 



n—l \ i—1 



Then, one easily writes 

Cit) = t'-i e-* • - I -V{St >y)~ x-^-^dx 

+ e-U-^{V{St + ii>y)^nii>y)} + 0{t^-^). 

Using the self- similarity of S, the second term on the right hand side converges 
to ip{y~) - p{y^)) E and so, for 1 < a < 2, 

limC(t) = {piy-)-Piy^)) ^ 0. 
4. Expansions via approximations by smooth functions 



The identity (1.6 1 suggests the possibihty of achieving power expansions for 
V {Xt > y) by approximating f(x) — l{3.>j,} using functions in Cf'+^ To 
this end, let us introduce moUifiers ipm G C°° with compact support contained 
in [—1, 1] that converges to the Dirac delta function in the space of Schwartz dis- 
tribution. For concreteness, we take iprn{x) := rrnp(mx), where 1^9 is a symmetric 
bump like function integrable to 1. Notice that 

ipm{u)du, (4.1) 

00 



converges to /(x), for any x ^ y. Clearly, applying (1.6 1 to each fm, 

j-k fn-\-l rl 

E/„(Xt)- VTTi'=/™(0) + — / (l-a)"E {L"+V,„(X„0}rfa, (4.2) 
k=i- ^" 

and by the dominated convergence theorem, 

lim EfmiXt)^V{Xt>y). (4.3) 

771 — ^OO 

Thus, the problem is to identify conditions for the limit of each term on the 
right-hand side to converge as m ^ cxd and to identify the corresponding limiting 
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value. The advantage of working with (4.2 1 instead of the decomposition of X 
of the previous section is that the coefficients dk of (1.3 1 can be identified more 
explicitly. 

As before, Cg £ denotes a smooth truncation function such that 
l[-e/2,£/2](a;) < c,{x) < l[„j^^](a;). 
The following operators will be useful in the sequel 

Lsgix) ■■= [ g{x + u)ce{u)v{du) 

g"{x + /3w)(l - I3)dl3w'^c^{w)v{dw), 



Lig{x) j 



where Ce{u) := 1 - c^{u), := -/ Ce{u)v{du), bi := 6 - / m(1{|„|<i} 
Ce{u))v{du) and 62 := Note that 

5 

Lg = ^ L,g, 



for any bounded g E C^. Moreover, it turns out that the following commuting 

LiLjg = LjLig. 



b ■ 

properties hold true for any g £ C^. 



Remark 4.1. Under additional assumptions on the Levy triplet {cr^,b,v), we 
can choose more parsimonious decompositions of the infinitesimal generator. For 
instance, if one of the hi 's is zero, then the corresponding operator is superfluous 
and can be omitted in the analysis below. Also, if Ji^i^i < +00 (i 



m 



which case the Levy process has bounded variation), then L/^ can be defined as: 
^4.9-= / {g{x + w) — g{x)) c^{w)v{dw) = I I g' {x + j3w)d(iwCi.{w)v{dw)^ 



provided that bi is adjusted accordingly. If t^(M\{0}) < 00, L4 can be omitted, 
provided that we define ^3,60 a^c^ ^1 via: L^g{x) — J g{x + u)iy{du), bo = 
r/(K\{0}) and 61 = 6 — uv{du). 

Let us introduce some more notation. For k := (fco, . . . , ^4) with fco, . . . , ^4 > 
0, u := (ui,...,Ufe3), w := (wi, . . . , WfcJ, and /3 := (^1, . . . , /S^g), define the 
finite measure 

d7r^(u, w,/3) = ]^Ce(Mi)i/(dui) J]^ c^{wj)wjiy{dwj) ]^(1 - f3j)df3j, 
1=1 j=i j=i 
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on the space -Ek ■— M'''3+''* x [0, l]*^*. Consider also the following related finite 
measure 

fcs ^4 ^4 

d7r^(u2, . . . ,Ufc3, w,/3) = ]^Ce(ui)i^((iuj) J]^ Ce{w-i)w^i^{dwj) ]^(1 - Pj)dPj. 

1=2 j = l j=l 

We sometimes drop the subscript k and superscript e in the measures defined 
above. Also, the integral of a function g with respect to a measure tt^ is denoted 
by 7r^(.g) and we assume, by convention, that T^^ig) = g, when ^ k4 — 0. 
Similarly, T^^ig) = g if = and fcs = 1 or 0. 

Let ICk be the class of all k = (fco, . . . , ^4) is such that fc^ > and fco+- • •+fc4 = 
k. Note that, for any fc > 1, 



(4.4) 



where (j^) — fc!/(fco! . . . k^l) is the multinomial coefficient and 



B^g{x) I ^C^.+^fe.+zM L + + ^/J^.^J dn, 



=1 3=1 



We first show that all terms in the right-hand side of (4.2 1 converges. 

Proposition 4.2. Let y > 0, n > 1, and < e < j//(n+ 1) A 1. Assume that v 
has a density s such that for any < k < 2n + I and any S > 0, 



ak,S '■— sup |s*-'^-'(a:;)| < 00. 

\x\>S 



(4.5) 



Then, for any 1 < k < n, 



rffe(y) := lim LVm(O) = V Ck(, ]a^, 



(4.6) 



where, for k — (fco, . . . , k^) and i?k := + 2fc2 + 2/24, Ck anrf Ck := ak(j/) are 



Ck 



ttk 



/(c,s)(^''-i) I 2/ - - I d%^, ks > 0, 4 > 0, 

i=2 j=l 



(4.7) 



A:3> 0,4 = 0, 
otherwise. 
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In particular, the limit 



TZn{t,y):= lim / (1 - a)" E L"+i/m(^at)rfa, (4.: 



exists, and moreover. 



n+l 



fe=l 



t" t 

fc! ^ ^ 



(4.9) 



Proof. We write dk for dk(y) and TZn{t) for 7?.„(t, y). From (4.4), it suffices to 
show that for any k = (fco, . . . , ^4) € /C^, 



lim Bkfrnix) = (-1) 



(/ci-l)l{f>0} 



Ok, 



with := A:i + 2fc2 + 2fc4 (recalling that by convention /.gdTTk = 5, if k^ + k^ = 0). 
In case £ — 0, 



which clearly converges to J '^{'^i^i^i > y} dn^, since fmix) l2:>j; for 

any x ^ y, \fm\ < 1, and TTk is a non- atomic finite measure. 

Consider the case k^, (. > Q. Writing z = X]i'=2 Pj'^j ^'^'^ integrating 

by parts, 

fm {x + ui + z) Ce • s{ui)dui = / y'm^^^ {x + Ui + z - y) ■ s{ui)dui 



= (-l)'^i~i y ^^{x + ui + z-y) {c,s)'^'^-^\ui)dui 
^-^{^lf^-\c,sr-^\y^x^z), 
provided that c e C^^^(M\{0}). Moreover, we have that 

fm {x + ui + z)c- s{ui)dui 



< 2^ ^ max sup |s*''''(a;)| < 00. 

k<£~l 

Thus, applying first Fubini's theorem and then the dominated convergence the- 
orem give: 

(c,s)(^-i)(y - X - V - V P,w,)d^ 

A O A 1 



i=2 i=l 



In case k-i = Q and £ > 0, 
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for m large enough since, by construction, e is chosen so that y — (n + l)e > 0, 
and (3jWj takes values in [—£,£] on the support of tt^. Then, the existence of 
and the identity (|4j| follow from (lO) and (|4!3|. □ 



Notice that we cannot yet conclude that dk :— d\^{y) are the same constants 



as the dfcS in equations (3.101 and (3.111, since we have not shown that 

limsup |7?.„(i, < oo. 



(4.10) 



Actually, in view of Theorem |3.2[ the se two conditions are equivalent; namely, 
dk = dk for all fc < n if and only if (4.101 holds. We show these facts in the 
following result. 



Theorem 4.3. Unde r the conditions of Theorem 3.2, (4^.10) holds and more- 
over, dk — dk in ( 4-6), which are independent of e and given by 

ri-l 



lim 
t— t 



1 



'{Xt>y)-J2^' 



fe=i 



fc! 



dn 

nl ' 



for k < n. 



Proof. Using a proof as in Theorem 3.2 we conclude that 
EfUXt) = e-^^* ^c,-„ - + Oe(i"+'), 

3 = 1 



(4.11) 



■j,m — ELj(i)^i 7i,m(0), with f^^„l{x) / /,„(x + u)s**(w)du. As 



where e 

Remark 3.3 (i), (|4.1l| leads to 



k=l 



with 



^k,m 



(4.12) 



(4.13) 



Since the last term in (4.2 1 is Oe(i"+^), we have 

dk,m = fm{^)- 



To show that dk '.= lim. 



,dk.m is identical to dk of (3.10)-(3.11), it suffices 



that lim,,„^oo Cj,m = Cj, or equivalently, that 

lim L,^/v^(0) = Lj/.(0), 



for all fc > and i>l. The case fc = is clear. For fc > 1, from (3.9 1, we only 
need to have 

Jhn^ j f^% L + Y.P^wA dnl = I f^ L + Y.PiwA dnl (4.14) 
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for any fc > and p > 1. Since ipm is symmetric, 



fi.m (•^) 



ipm * S*\u)du, 



and thus, ,„(a;) = ipm * sl^iu — 2;), which converges to //(a;) = — x), as 
m — > 00, uniformly in x. Then, 



which converges to 



{-l)P-^st^^-^) , 4f-i)(y -x)^ fP\x), 



as TO — !■ 00, uniformly in x. Since 7r| is a finite measure, (4.14) holds. We have 
just proved that dk = dk, for k < n, and by matching (3.101 and (4.9), it follows 
that TZn{t) = 0(1), as t 0. The last two statements of the result are easily 
proved by induction. □ 



5. The remainder and expansions for the transition densities 



In this section we give a more explicit expression for the remainder TZn in (4 



whose existence was proved in Theorem |4.2| In order to do this, we expand 
L^~^^ fm{Xat) using (4.4) and show that the limit of the resulting terms exists. 
The hardest case to tackle corresponds to — and £ > 0, where we will need 
to impose the following condition on the transition density pt of Xt : 



Cfc A 



sup sup \pI^\x) \ < 00, 

0<«<to \x\>S 



(5.1) 



for any S > and for some > 0- Condition (5.1) is reasonable since it is 
known that 

-ptix) - s{x) 



lim sup 

"|a;|>£ 



t 



0, 



(5.2) 



(see Proposition III. 6 in |T^ and Corollary 1.1 in [HI)- We confess however that, 
in general, this condition might be hard to verify since the transition densities 
Pt of a Levy model are not explicitly given in many cases. Let us point out that, 
under certain conditions, Picard |14j proves that 



sup \p^^\x)\ 



< t 



-(k+l)/0 



(5.3) 



where (3 is the Blumenthal-Getoor index of X. The approach in [14 was built 



on earlier methods and results of Leandre [12], who proves (5.2) and (5.3) for 
A; = using Malliavin calculus. In view of ( |5.3p , for values of t away from 0, 
the derivatives of pt are uniformly bounded, and condition (|5.1| is then related 



to the behavior of pi'"' when t ^ 0. In Sections [6]and|7l we prove that the 
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condition (5.11 holds for symmetric stable Levy processes and other related 
processes, rising the hope to use similar methods in other cases. As in the 
previous section, we take y > 0, n > 1, and e > such that 



< e < y/{n+ 1) A 1. 



(5.4) 



Theorem 5.1. Assume that v has a density s such that {4-5) holds for any 
< fc < 2n + 1 and any d > 0. Also, assume that there exists a to > such 
that for all < t < to, Xt has a (7^"+^ density pt satisfying {5.1) for any 
< fc < 2n + 1 and any (5 > 0. Then, the remainder 

TZn{t,y):^ lim / (1 - a)" E L"+ V^(X„t)(ia, 



is given by 



(1 — a)"ak(i; a, y)da, 



(5.5) 



kG/C„+i 

where, for k = (fco, . . . , k^) G K-n+i, Ck, and ak{t; a) are defined via: 
Ck:=6^''6t^6^H-l)(^-i-i)i{^>">, 



JP[X^t + J2u,>y\dn^ 



£ = 



i=l 



ai^{t;a,y) := < 



y-X^t- Y,^, - Y,l3,wj \d%^, fcg > 0, £ > 
/pit'^ I y - Y^Pj^i I dn^, fcs = 0,£ > 0, 



with £ := ki+ 2fc2 + 2fc4 



Proof. From (4.4 1, it sufhces to show that for any k — (fco, . . . , k^) € ICk, 

lim /V-«)"Ei?k/„.(^at)rfa=(-l)''^"'^'^'>°^ [\l~ara^{t;a)da. 

Below, Ti^g{x; •) is the function defined on Ei^ := M'^^+'^i x [0, 1]'''' via 
Tkg{x; ui,. . . ,Uk3,wi, . . . ,Wki,Pi, . . . ,/3feJ 



^ik^+2k,+2k,) L + + 



We break our proof in different cases. Suppose first that £ := fci +2fc2 + 2fc4 
0. Since < /m < 1, apply Fubini's theorem to get: 



EBMXct) = 6g"6^^6^^^^(ETk/,„(X,t, •)). 
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Since X^t is a continuous random variable, the dominated convergence theorem 
impUes that 



2 = 1 



Again, by dominated convergence, /g'^(l — aY^ E B]g_f„-i^{Xat)d<^ converges to 
/o (1 - a)" / P >y- u,) dTT da. 

Next, we consider the case i > and = 0. Again by Fubini's theorem. 

Writing z = X]j=i Pj'^j^ integrating by parts, and changing variables, we have 
Efg^ (X„t + ^)= J ^fn^^ (a; + z - y)pat{x)dx 

which converges to {—\)^'^~^p'^^ — z) as m — s- oo, if p^^^ is continuous. 



Moreover, under (5.1) and with the help of (5.4 1, for m large enough 



sup sup 

0<a<l 0j,Wj 



Ef^UX^t + Y,P, 



< sup sup Ip^^t ^\x) \ < oo, 

0<a<l \x\>6 



taking S := {y — {n + l)e)/2. Then, by dominated convergence: 
lim / (l-a)"EBk/™(Xat)da = 

Note that the previous limiting value is uniformly bounded in t and y by 
^7r(M2'=4) gup sup < oo. 



n + 1 



0<u<t \x\>S 



The only remaining case to tackle is when £ > and > 0. Writing z = 
Si^=2 ^« + Si=i have that 

/i^^ (^at + ui + z) {cs){ui)dui = y ^(,^^^1) {X^t + ui + z-y) {cs){ui)dui 

- f ^r,AXat + Ui + Z - y) {csY''^\ui)duu 
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which converges to (—1)*^^ ^{cs)^^ — X^t — z) as m ^ oo, provided that 
c e C^^^(M\{0}). Moreover, under (4.51, we have that 



f^^ {Xo,t + + z) {cs){ui)dui 



< 2^ ^ max sup \s''''^x)\ < oo. 

k<i-l \x\>e 



Thus, applying first Fubini's theorem and then the dominated convergence the- 
orem give: 

lim / (l~a)"EBk/™(X„t)da = 



1=2 



This last case achieves the proof of the theorem. 
Remark 5.2. From the proof is clear that 



□ 



\T^nit,y)\ < 



an 



n+1 



max < max Cfe.^, max a^g,! 



k<2n+l 



k<2n+l 



for any t G (0,to); where 



an 



ke/c„+i 



n+1 
k 



max{7rjl),^i-jl),7r^(l)}. 



This hound on Tl{t, y) is valid for any y > y, taking also e such that < e < 
y/(n+l)Al. 

One of the advantages of an explicit expression for the remainder TZn{t, y) of 



(4.9 1 is that we can obtain small time expansions in t for the transition density 
Pt{y) of Xt, by a formal differentiation of ( |4.9| . With this application in mind 
we need to show that the coefficients ak(y) := ak of (4.7 1 and 



^k(t,y) 



(1 - a)"ak(i; a,y)da, 



of (5.5 1 are differentiable in y. The following result corrects Theorem 1 from 

Proposition 5.3. Let y > 0, n > 1, and < e < y/{n + 1) A 1. Assume that 
V has a density s such that (4--5) holds for any Q < k < 2n -\- 2 and any 5 > Q. 
Also, assume that there exists a t^ > Q such that for all < t < tQ, Xt has a 
(j2n+2 dg^gHy satisfying (5.1) for any Q <k <2n + 2 and any (5 > 0. Then, 

n J. 

My)^y^dUy)Ti + Oe{r+'), (5.6) 



k=l 
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where for any I < k < n, 



d'kiy) ■■= Ck(, )ak(y) 



(5.7) 



and, for k — (fco, . . . , k^) and £k '■= ki + 2^2 + 2fc4, 



(^3 ^4 \ 

y - - E^j^'j ^3 > 0' 4 > 0, 

i=2 3=1 I 

/ fc3 \ (5.8) 

J{ces) [y -^uA dn^, A:3> 0,4 = 0, 



0, 



othe 



Proof. In view of (4.9 1 and (5.5), we first have to show that the derivative 
of (4.7 1, for each case, exists and is given by (5.8 1. This will follow from the 
fact that TTk and TTk are finite measures and the integrands have continuous 
uniformly bounded derivatives with respect to y. Also, we need to show that 
the derivatives, with respect to y, of ak(t; a, y) exist and are continuous as well 
as bounded for t small. When £ = and fca = 0, 



Ai,{t,y)= (l-a)"ak(i;a,?;)= / (1 - a)" / Po,t{z)dzda. 



From (5.2), there exist K > and to > such that supq^^^j^^ sup^^^ygPu{x) < 
K. Hence, one can interchange derivation and integration: 



dA],{t,y) 
dy 



(1 - a)'^pat{y)da, 



(5.9) 



and moreover, the supremum of (5.9) over Q < t < to, is finite. In case of ^ = 
and /ca > 0, one can write 

^k(i, j i^-'^T j J Pat \ z - Yj^'J dT^kdzda. 
The inner integral is continuous in z and is such that 

< Patiz) := J ■ ■■ J Pat ^z - YJ^^ {CeS){ui)duidTTk < supSe(ui)A^'='~^ 

Thus, one can interchange derivation and integration to get: 
5Ak(i,y) 



dy 



(1 - a)'^pat{y)da, 
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and moreover, the suprcmum over t > is finite. In case > and £ > 
fi 



^3 ^4 

(1 - a)" / E (cs)(^-i) [y-X^t- - Y^f3jWj I dn^da. 



i=2 



Assuming that cs € C^, one can interchange the derivative with respect to y 
and the integral, and the resulting term will be uniformly bounded in t since TTk 
is a finite measure. Finally, when fca = and ^ > 0, 



y) / V - «)" / pi'r'' iv - E^.^. j dn. 



da. 



Assuming that p^t G and satisfies (5.11 with k = £, one can interchange 
the derivative with respect to y and the integral, and the resulting term will be 
uniformly bounded in < t < to since TTk is a finite measure. All previous cases 
will imply that 

dAi,{t,y) 
dy 

exists an d is uniformly bounded in < t < to- Hence, the derivative of the last 
term in (|49| is 0^{r+^). □ 



6. Symmetric stable Levy processes 



In this section, we analyze the assumption (5.1), needed for the validity of 



Theorem 5.1 and |5.3| in the case of symmetric stable Levy processes. 

Let us assume that the Levy triplet (a^, b, v) is such that 6 = and that v 
symmetric. Furthermore, let us assume that 



Km inf 



> 0, 



(6.1) 



for < a < 2. Condition (6.1 1 is equivalent to 



lim e° 



{\x\>e} 



v{dx) > 0. 



Condition (6.1 1 is known to be sufficient for Xt to have a C°°-density pt (see 



e.g. [m Theorem I.l] or [19, Proposition 28.3]). It will be useful to outline 
the proof of this result. The first step is to bound the characteristic function 
iptiu) = Ee™^* as follows: 



|^t(w)|<e-^*H°, 



(6.2) 



which is valid for u large enough (cf. page 190 in [19 ). Note that the right hand 



side of (6.2) is the characteristic function of a symmetric a-stable Levy process. 



In particular, 



/ 



|V't(u)| \u\'^du < oo, 
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for any n — 0, . . . , and the following inversion formula for pj"-* holds: 

= ^ / e-^-VM^)du, (6.3) 



27r 

see |19l Proposition 2.5]. Finally, the Riemann-Lebesgue lemma implies that 

lim p["\x) — 0. 

\x\ — »oo 

Let us try to modify the above argument for our purposes. In the case that 
6=0 and i' is symmetric, iptiu) is positive real and even, and thus, 

- — ^ — / cos{ux)u"'ipt{u)du, if n is even, 
pi''\x)^\ ° (6.4) 

sm{ux)u"'ipt{u)du, if n is odd. 



(_l-)(n + l)/2 



In light of (6.2 1, it is important to analyze the case of a symmetric a— stable 
Levy process. It is not surprising that a great deal is known for this class (see e.g. 

Section 14 in [IE,). For instance, from the self-similarity property Xt ^t^^^Xi, 



Pt 



Asymptotic power series in x are available for pi(x), from which one can also 
obtain the following asymptotic behavior of pi{x) when a: ^ oo: 

Pi(:e)~x-"-\ (6.5) 

Note that (6.5 1 is consistent with the well-known asymptotic result that 

1 



lim - Ptix) — six) = x' 
for any x (see e.g. [13 Corollary 1]). 



-Q-l 



We want to show that the condition (5.1) holds for symmetric stable distri- 
butions (and possibly for more general symmetric distributions satisfying (6.1 1). 
With this goal in mind, we give a method to bound x'^^'^pi{x). First, we need 
the following lemma: 

Lemma 6.1. Let (f> : (0, oo) M_(_ be an integrable Junction. Then, the following 
statements hold: 

(i) If (f) is monotone decreasing and there exists (3 G [0, 1] such that 

limsup 0(m)m'' < oo, (6.6) 



then there exists a constant c < oo, independent of x, such that 

K{ux)(f){u)du 
where n can be either cos or sin. 







< (6-7) 
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(ii) If (j) is unimodal with mode u* , then 

K{ux)(j){u)du 



< 



Tr(f>{u*) 



(6.8) 



for all X > 0, where k can be either the function cos or sin. Moreover, if 
(j) is continuous, then 



(6.9) 



lim X / K{ux)(j){u)du — 'K(j){u*). 
Jo 

Proof. To show (i), we first note the following two easy inequalities: 
< / sm{ux)(i>{u)du < / (j){u)du < oo, 



cos{ux)(f){u)du 



< 



(u)du < oo, 



(6.10) 



(6.11) 



valid for any nonnegative function (f> that is decreasing and integrable. Therefore, 
if K, is either cos or sin, then 



ti{ux)<j){u)du 



< 



4>{u)du. 



In view of the condition (6.6 1, there exists a c' > and xq > su ch that for 
all X > Xq, (j){u) < c'u"!^, in (Q,2n/x\, for any x > Xq. Then, (6.7) is clear for 
X > xq. The values x < xq can be taken care of easily since 

1^2-17 / X poo 

(p{u)du t / (j){u)du, 



when a; \ 0. Let us now show (ii). First, set 



q{x) 



K(ux)(j){u)du. 



By assumption, (j) is increasing on [0,w*], and decreasing on [m*,cx)). It can be 
shown that for any x > 0, there exists a positive number u{x) such that 



k{xu{x)) = 0, \u* — u{x)\ < 



2n 



u(x) 



u{x)~-K I X 



n{ux)(j){u)du < q{x) < 



and 

'u(a;)+7r /x 



u(x) 



K{ux)(j){u)du^ 



(6.12) 



(see e.g. Figure [T] where the choice of u{x) is illustrated when k{u) = cos(w)). 
Next, the upper and lower bounds on q are such that: 



(x) 

u{x) — 7t/x 



ti{ux)(j){u)du < — (l){u{x)) < —(j){u*), 

X X 



TT TT 

tl{ux)(j){u)du > (j^iuix)) ^ 
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where u{x) € [u{x),u{x) + tt/x] and u{x) G [u{x) — Tr/x,u{x)]. The inequality 
(6.81 is thus clear, while (6.9 1 results from the fact that both u{x) and u{x) 
converges to u* as a; ^ 00. □ 



25 



'□ 10 20 30 40 50 60 



Fig 1. Definition of a{x) 



Proposition 6.2. Let X be a symmetric a-stable Levy process, and let pt be 
the density of the marginal Xf . The following two statements hold: 

(a) //O < a < 1; then there exists an absolute constant c such that 

sup \x\'^^^pi{x) < c. 

X 

(b) Lf I < a < 2, then for any e > 0, there exists a constant < c(e) < 00 
such that 

sup |a;|"+V(2;) < c(e). 

\x\>e 



Proof. Without loss of generality suppose that x > 0. By (6.4 1, the well-known 
representation of the characteristic function of Xt, and an integration by parts, 

Pi{x) = — cos{ux)e~^ du — — sm{ux)u°'~^ du. (6.13) 



TT 







X 







If < a < 1, then we can apply (6.7 1 with (3 = \ — a, and hence. 
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for a constant c. Now, let 1 < a < 2. Applying another integration by parts in 

a{a — 1) 



(6.13 1, we have 



cos{ux)u°' " du 
cos(ux)M2("-i)e-"°dw. 



(6.14) 
(6.15) 



The first term in the previous inequality can be bounded using (6.7 1 with /3 
2~a: 



1 



cos(mx)u° 



du 



< 



,1-/3 



ra+1 ■ 



The term in (6.151 can be bounded using (6.8 1 since (j)(u) = u 



2(a-l) 



IS 



unimodal and thus 



< c-T < c'— — . 

— X — X ' 



for all X > e, where c, c' < oo are constants depending only on e. Plugging in the 
above bounds in (6.151, we obtain the second statement in the proposition. □ 

Remark 6.3. In view of the above proposition, we obtain the following bound 
for the transition density pt of a symmetric a— stable Levy process: 



pt{x) < c- 



t 



2,Q + 1 ' 

valid for all t > and \x\ > e, and where c is a constant depending only on e. 



We can now generalize the ideas of Proposition 6.2 to dealt with the deriva- 
tives of the transition density. 



Theorem 6.4. Under the conditions of Proposition \ 6.2\ for any e > 0, there 
exists a constant c„(e) such that 



sup \x 

|a:|>e 



p[''\x) <c„(e). 



Proof. We prove the following more general bound: 



sup \x 

\x\>e 



a+l+n 



K{ux)u"-e du 



< dn{e) < oo, 



(6.16) 



(6.17) 



where k can be either cos or sin. Without loss of generality, let us assume that 
X > 0. Our proof is then performed by induction on n. Proposition |6.2| yields 



(6.171 for n = and k{x) = cos(a;). The case k{x) = sin(a:) can be dealt with 
in an analogous way; namely, we first integrate by parts, once when a < 1, or 
twice when 1 < a < 2, and secondly, we use (6.7 1 if a < 1, or (6.8) if 1 < a < 2. 



J.E. Figueroa- Lopez and C. Houdre/Small-time expansions for Levy processes 



29 



Now, assume that (6.171 holds for n = 0, . . . , m — 1. We want to prove the 
case n = m > 1. Set 



qni{x) 



K.{ux)u™'e "^"du. 



Applying consecutive integrations by parts, one can find constants bj (depending 
only on a and m) such that 



1 



K(ua;)w"~^e~" du + 



— y"^j / K(wa;)u*"e-"°(iM, (6.18) 
^" ,=1 Jo 



(6.18 1 is such that 



where k, k are either cos or sin. By the induction hypothesis, the first term in 

K(Mx)u™-ie-"°(iM 



sup 

\x\>e 



1 

xJo 



(6.19) 



as we wanted to show. 

Now, for the second term, let us consider first a < 1. Let fc > 1 be such that 

fc- 1 k 

< a < 



k — 1 + m k + m 

Also, for each 1 < j < m + fc, let 1 < < j be such that 



r, - 1 



< a < 



3 



J 



Setting S{x) := X^JLi Io° i'^{ux)u^°'e du, and applying successive integra- 
tions by parts to each of the terms of S{x), it follows that 



m+k „oc- 

.7=1 



du 



(6.20) 



for some constants a^, and where Kj is either cos or sin. By the way Vj is chosen, 
the inequality (6.7 1 can be applied to estimate the absolute value of each term 
in (|6.20|. Then, 



K,j{ux)v?°' ''^e " du 



(6.21) 



for some dj > 0. Combining (6.18l-(6.21 1, there exists a constant Cm(e) such 
that 



sup 

|a;|>£ 



p(")(x) |xr+l+'"<C„(£). 



(6.22) 
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Next, we consider the case of 1 < a < 2. Note that, for some constants oq, Oi, 02 
depending only on a and j, the term Cj{x) := Jp K(ux)u^°'e^" du can be 
broken into three pieces: 



Ol 



k{ux)u 



02 



(6.23) 



If J = 1, then the first term in (|6.23| can be bounded using (6.7) with (3 = 2 — a: 

r. r 

< 



k{ux)u 



ra+1 ' 



for some c < 00. For the other two terms of the case j — I or any other 
2 < i < Ti, we can apply ( 6.8 ) since then the function multiplying k is unimodal. 
Then, for any e > 0, we can bound <S'(a;) :— J2'f=i bj Jq K,{ux)u^"e^'^ du in the 
following way: 



sup|S'(x)| < 



c c 
— < 



for a constant c" depending only on s. 

Finally, let us verify the case a — 1. Without loss of generality, assume that 
k{x) = cos(a;). After two integrations by parts, we have that 



cos{ux)u"^e ^du 



x 
m 

X 



sin(ua;)u™ + 



sin(ua;)u" 



'du- 



sin(ux)u™e "du 
cos{ux)u"^ du 



+ ^ cos(ua;)u"~^e^"du 



We can then write the above equality in the following manner: 



1 



cos(ua;)u™e "du ■ 



sm{ux)u"^ "^du 



m 



cos{ux)u 



The result follows by applying our induction hypothesis to bound each of the 
two terms in the right-hand side of the last equality. □ 

Corollary 6.5. With the notation of Propositioned^ for any < a < 2, e > 0, 

and n > 0, there exist a constant Cn.e such that 



sup 

\x\>e 



pi-\x) 



(6.24) 



for any <t < 1. 
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7. General Levy processes 



In this part, we examine the vahdity of the assumption (5.11 for general Levy 
processes, whose Levy density s is stable like around the origin. 

The main tool will be a recursive relations between the derivatives of a density 
p. Consider a distribution /i such that its characteristic function ip{u) := fl{u) 
is C°° with also 

/•OO 

iP^^'^u) du < OO, (7.1) 

for all r > and m > 0. Recall that in that case fj, admits a C°°-density p and 
moreover, 

pM(2.) ^ ir^r f e-™-u™^(u)dM. (7.2) 
2n J 

By applying two consecutive integration by parts, we can derive the following 
formulas 

where we are assuming that m > 2. However, even if m < 2, we can deduce 
a recursive formula for terms of all its lower order derivatives and the 

integral of the function 

du<' 

Indeed, we have: 

Theorem 7.1. Let r > and m > 0. Then, for all x, p^"^\x) can be written 
as 

YC, tf(m- z)lp("-^-)(x) + {^lytTZL L-^u.^,nd[jMdu, 
where c™- are given by the following recursive formulas: 

Proof. We prove the formula by induction in m. Consider the case m = 0. We 
want to prove that 

.(x) = (-l)^(^/e— |,^(.)d., 

for any r > 0. This can be done by induction on r and integration by parts. 
Suppose that the formula is valid for m — k and all r > 0. We want to show the 
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formula for m = fc + 1 and all r > 0. Now, we use induction on r. The case r = 
is just (7.2 1 with m = k + 1. Suppose the result holds for r = £ and m = k + 1: 

i/\{k+l) 



J-1 



P 



(k- 



3=1 4=0 

\k + l-l 



P 



(7.4) 



Next, with an integration by parts in the last term 

M(fe+1) j-1 



j = l i=0 

\k+i-e-i 



\e+i ( 



(7.5) 



Then, writing (7.3 1 for m — k and r ~ £ and solving for the last term gives 



(-1)^ 



k-e 



du^ 



(^) + E<. T[{k-^)^P^'-'\x) 



Plugging in (7.5 1, we get (7.3 1 with r — I + \ and m = fc + 1 provided that we 
define the coefficients c^+i ^ as follows: 



I I - _ 1 ^fi^+J- ^K+i I K 



This proves the case of r = ^ + 1 and so, the result holds for all r and all m. □ 

The following corollary give further information when working with the tran- 
sition distributions of a Levy process. 

Corollary 7. 2. L et (Xt)t>o be a Levy process such that fi, the distribution of 
Xi, satisfies (7.1). Let 7 be such that i^ti^) :— e*'^^^\ where tpt is the charac- 
teristic Junction of Xf. Then, the density pt of Xf admits the representation: 



E 

3 = 1 



.{-^r 



p'r\-) = E <3 IK- - ^)i^p'rH-) + i-iy 2,,. 



■^"\t,x), (7.6) 



i=0 



where 



/or some constants d^^'j^ . The above summation is over all non-negative integers 
*iji2, ji,i2 such that < ^2 < ii and iiji + 12^2 — r. 
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As an application let us consider a Levy process as in Corollary |7.2| such that 
for each « > 1, there exists Ci < oo and uo,i > such that 

M*H^)l<c.H"-\ (7.7) 
for all \u\ > uo,i. Also, assume that there exists Uq > and cg < oo such that 

<e-^°l"l°, (7.8) 



for all u > uq. Remember that (6.1) implies the above condition (cf. Sato [HI 
Proposition 28.3]). Then, we have the following result: 



Proposition 7.3. Let (7.7) and (7.8) be true for < a < 2. Then, for any 
m > 0, any e > 0, and any tQ > 0, 

sup sup |p(™^(a;)| < oo. 

0<t<*o |x|>E 



Proof. The proof is by induction on m > 0. The recursive formula (7.6) with 
r = 1 and m = leads to |pt(a;)| < t \J e^^"'^y{u)e*''^'^")du\ jx. Note that we can 
assume that there exist constants mq > 0, 6o, and h\ such that 

sup |7'(w)||e*^('")| < 6o, \i{u) ■ e*T(")| < 6i|u|"-ie-=«*l"l° , 

|l«|<«0 

for all \u\ > Uq and < t < tg. Then, for all t < tQ, 

t r°° 

\pt{x)\<bQUQ- + b, / v'^-^e-^°^ dv. 

X Jq 

Next, let the statement of the proposition hold true for m = 0, . . . ,k, and let us 
show it for m = k + 1. In view of (7.6), it suffices to show that 

sup sup |I™(t, a;)| < oo, 

0<t<to |£c|>e 

for some r > 0. Moreover, it suffices to show that 



sup sup t^^'^^^ 

0<t<to \x\>e 



< OO, 



for any ii > ^2 > and ji, j2 > such that iiji + 12^2 = f. As before, we can 
assume that there exist constants mq > 0, &01 smd bi such that 

sup |7('i)(M)P'i|7('^)(u)P'^|e*''(")| < 60 

|u|<tio 

for all \u\ > Uq. We need to show that there exists an r such that the supremum 
on < t < io of 

/•oo poo 
t'h+h I y(il+j2)Q+m-rg-cotu°^^ _ ^i(r-m-l) / y(jl+h)a+m-r ^-cqv" 

Jo Jo 
is finite. The supremum above will be finite if r = m + 1. □ 
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Example 7.4. Consider the CGMY Levy model introduced in 121 and of great 
popularity in the area of mathematical finance. This process is a tempered stable 
one in the sense of Rosinski llSf . Its characteristic function is given by 

Mu) = exp{tCT{-a) ((M - iu)°' - M" + (G + m)" - G")} 

(see Theorem 1 in 121). Then, 

7(u) -.^ CT{-a) {{M - iu)" - Af" + (G + m)" - G") . 



We can then verify that 7 satisfies (7.1) and ( 7.( 



The next result generalizes the conclusions in the above example to more 
general tempered stable processes. For simplicity, we take symmetric processes, 
even though the proof can be extended to the general case. 

Proposition 7.5. Let X be a Levy process with Levy triplet (0,0,;^). Assume 
that V is of the form v{ds) = |s|^"^^(7(|s|)(is, where < a < 2 and q is a 
completely monotone function on M-|_ such that 



q{s)ds < 00, 



(7.9) 



for all j > I. Assume also that the measure F for which q{s) — e ^''F{dX) 
is such that 



X^F{dX) < 00, 



(7.10) 



for all j > 0. Then, the function 7 associated with the c harac teristic function of 
X via i^tix) '.= e*''''"-' satisfies the conditions (1.1) and [1.8). 



Proof. Clearly, 



-2-a 



(7.11) 



and thus, condition (7.8 1 will follow. Now, we claim that there exists a constant 
C such that 

/•OO 

< Gu"-i 



sin(us)s 



(7.12) 



for all A, w > and < a < 2. Indeed, if < a < 1, (|712| results from (|6/7f . If 
1 < a < 2, then changing variables and using sinz; < v, 



sin(us)s "e ^^ds 



sm{v)v-"e-^'"/''dv 



v^-°'dv + u°'-^ 



y-^dv < Cu° 



for a constant C independent of u and A. Moreover, it can be proved that there 
exists a constant d such that 



K{us)s^-°'e-^''ds 



(7.13) 
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for J > 1, A, M > 0, and < a < 2, and where k can be either cos or sin. Indeed, 
the case j — I can be proved as fohows. If < a < 1, then we apply two times 
integration by parts (similar to the case a = 1 in the proof of Theorem 4.4). 
Then, we can apply part (i) of Lemma 6.1 If 1 < a < 2, then one can apply 
directly part (i) of Lemma 6.1 The case j > 1 can be proved using induction 
on j with the help of two integration by parts. From the previous estimates, we 
have that, for j > 0, 



j„.a-l-j 



(7.14) 



where Cj is a constant independent of A and u and Kj(u) — cos(ii) if j is odd, 
and Kj(u) — sin(w) if j is even. Next, from the conditions on X, the function 7 
is given by 7(1*) = 2 (1 — cosus)s^°'^^q{s)ds. Condition (7.9) implies that 



^0)(u) 



= 2 



Kj-i{us)s-' " ^q{s)ds 



J > 1, 



where above. In that case, using (^7.10|, applying Fubini's Theorem, and 



(7.141, we have 



< 2 



F{dX) 



< 2Cju"~^ / (1 + Xy-^F{dX) < Cju"-^, 



for a constant c, independent of u. 



□ 



Remark 7.6. Rosinski \16^ (see Proposition 2.7) gives conditions for {7.9) to 
hold. In terms of the notation of Proposition 7.5. (7.9) holds with j > I if 



and only if X ^ F{dX) < 00, which is also necessary and sufficient for j ~ 1 
provide d, tha t a < 1. If a > 1, then |y.ffp always hold for j — 1, while when 
a = 1, (7.9) hold with j = 1 if and only if X^^ log(A^^)i^(dA) < 00. 



Appendix A: Verification of the claim in Remark |3.3| (iii) 



Note that the expression for ^2 in Remark 3.3 (ii) can be modified so that one 
can replace Se{x) by s(a;)l{|2,|>gj.. We will get: 



d2 = -ah'iy) + 2bs{v) -2 f f s' 

J\w\<e JQ 



\x\>6 J\u\>E 



\w\<E . 

l{x+u>v}s{x)s{u)dxdu 



iy- I3w){l- (3)dl3w^s{w)dw (A.l) 



2 / s{x)dx I s{x)dx — 2 I xs{x)dxs{y). 

'\x\>e Jy ■Je<\x\<l 



(A.2) 
(A.3) 
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The last term in ( A.l I converges to as e ^ 0. The term in (A. 2 1 can be written 



as follows (omitting the integrand s{x)s{u) and using symmetry of this about 
X = u): 



z l*oo poc poo 

2 1 dx du + 2 dx I du- 

' — oo y~x J y 



OO /'OO 



dx / du 

y -'y 



A'2 

y ry ryl"^ ry 

dx I du + 2 I dx I du . 



y/2 •'y/2 



y-x 



As 



Similarly, we can decompose the terms in hne (A. 3 1 as 

-2 



s{x)dx I s{u)du — 2 I s{x)dx / s{u)du 

— 1 J y J — OO -/ y 



Bi 



B2 



-2 / s{x)dx / s(u)du—2s{y) / xs(x)dx —2s{y) j xs{x)dx 

ly Je J-1 



B, 



Bi 



B. 



Now, A2 + B3 = 0, Ai + Bi + B2 + Bi becomes 

/-e I'y ^-1 j'y 

I {s{u) — s{y)} dus{x)dx -\- 2 I I s{u)dus{x)dx, 
-1 J y—x J —OQ Jy~x 

and A5 + B^ becomes 

ry/2 py rl 

{s{u) — s{y)} dus{x)dx — 2s{y) / xs{x)dx. 

y-x Jy/2 

Then, after changing variables, d2 becomes: 

- cr^ s\y) + 2s{y)b + 2 J {s{y - u) ^ s{y)} dus{x)dx - y{[y, oo)f 

ry ry r-i ry ri 

+ / s{x)dx / s{u)du + 2 I / s{u)dus{x)dx — 2s{y) I xs{x)dx 

Jy/2 Jy/2 J -00 Jy-x Jy/2 

/-e I'X |.y/2 px 

/ {s{y ~ u) — s{y)} dus{x)dx + 2 / / {s{y — u) — s{y)} dus{x)dx. 
-1 Jo Je Jo 

Now, taking e ^ in the above and further manipulation, gives the expression 



in Remark 3.3 (iii) 
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